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Abstract 



The pion-nucleon coupling constant Qt^n is studied on the basis of the 
QCD sum rules. Both the Borel sum rules and the finite energy sum 
rules for g-,^N are used to examine the effects of higher dimensional oper- 
ators (up to dim. 7) and corrections in the operator product expan- 
sion. Agreement with the experimental number is reached only when 
S-k/Sn is greater than one, where S-j^ (Sn) is the continuum threshold 
for the g-j^N (nucleon) sum rule. 
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I. INTRODUCTION 



The pion-nucleon coupling constant g^^N is one of the most fundamental quan- 
tities in hadron physics. Qt^n = 13.4 ± 0.1 is obtained empirically using the N — N 
scattering data, n — N scattering data and the deuteron properties |]l[. From the 
theoretical point of view, it is a great challenge to reproduce this value from the first 
principle namely the quantum chromo dynamics (QCD). So far, there have been two 
attempts; one is based on the lattice QCD simulations (see ref. 0) and the other is 
based on the QCD sum rules (QSR) PJl]. The latter approach for g.^^, however, has 
not been explored in great detail beyond the leading order of the operator product 
expansion (OPE) within the authors' knowledge.^ The purpose of this paper 
is to reexamine the problem using the currently available information on the higher 
dimensional operators in OPE and the corrections to the Wilson coefficients. 

To study Qtjn in QSR, two methods have been proposed so far: (i) method based 
on the two point function (0 | Tr]{x)f]{0) \ vr), and (ii) method based on the three 
point function (0 | Trj{x)qij5q{y)f]{0) \ 0), where ri{x) is the nucleon interpolating 
field. We will take the first approach (i) throughout this paper, since OPE in (ii) 
is problematic when the four momentum of the pion becomes soft. The lowest non- 
trivial order of OPE in (i) is known to give g^°^'^^* = Mtv/Ztt when the continuum 
threshold is neglected 0,0. Th origin of the 25% disagreement of g^°^'^'^^ from the 
Goldberger-Treiman (GT) relation gf^^r = gA^N/ fir niay originate either from the 
higher dimensional operators, corrections and the continuum threshold. We will 
examine g-^N with all these ingredients. 

In section II, we will examine QSR for the nucleon mass by adopting OPE up 
to dimension 7 operators with 0{as) corrections, which is essential for the discussions 
in later sections. In section III, QSR for g^^^ is studied in close analogy with that 
for the nucleon mass. In section IV, Borel analyses for (^^at are made, and the effect 
of the Us corrections, higher dimensional operators and the continuum thresholds are 
studied. Section V is devoted to summary and concluding remarks. 

II. SUM RULES FOR THE NUCLEON 

Let's consider the two point function. 



There is a QSR study of the axial charge qa^ which can be related to g-j^N assuming the 
Goldberger-Treiman (GT) relation |^. In the present paper, however, we will study Qj^n 
directly. 
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na/3(g) = I J d^xe"^-^O\T7]^ix)f]0iQ)\Q) 

= U,{q)q^fS + U2iq)lafS, (1) 

with q = q ■ 'J. For interpolating nucleon current, we use the loffe current ^ 

V{x) = eaUu''{x)C^y{x))^,Yd%x), (2) 

where a, b and c are color indices and C is the charge conjugation operator. 
The operator product expansion (OPE) at — oo has a general form 

/ rf^xe*^-^T77(x)r/(0) = ^ /z, «,(/i2))0„(/i), (3) 

where C„ denote the Wilson coefficients and 0„ are the local gauge invariant opera- 
tors. They depend on the renormalization scale /i which separates the short distance 
dynamics in C„ and the long distance dynamics in 0„. If one takes the vacuum ex- 
pectation value of (0), it can be used for the sum rules of the nucleon, while if one 
takes the vacuum to pion matrix element, it can be used for the sum rules of Qt^^. 

The Lorentz structure of 0„ depends on the states one chooses to sandwich (|^). 
For the nucleon sum rules, only the scalar operators contribute. OPE up to dimension 
7 with as corrections in the chiral limit can be extracted from refs. [§] and ; 

Ren(Q^) = (n; + n; + n;) « + (n^ + n; + n^) , (4) 
n;(Q^) = ^onn9l + i^^im^i . (5) 

n?(Q') = T^(-G^)ln^, (6) 

n^(Q') = 0, (9) 

^kQ') = T^,{-G'){dd), (10) 

where = —q^ oo and (■) denotes the vacuum expectation value. The argument 
of as is /i^ which is not written explicitly. The diagrammatic illustration of 11" ~ li^ 
is shown in Fig. 1. 

Several remarks are in order here, 
(a) We take the chiral limit (mg = 0) throughout this paper. Small m, d quark mass 
does not change the essential conclusion of this paper. 
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(b) The fact that the Wilson coefficient of the dimension 5 operator gqa ■ Gq vanishes 
in 112 is a unique property of the loffe current 0. Since this operator is aheady 0{g), 
we do not consider the 0{as) correction to it. 

(c) A small discrepancy between the formula in [P] and that in |]T0[ for Ilg has been 
recently resolved (see [|lT|) and the final result boils down to the form given in the 
above. 

(d) We always assume the vacuum saturation when evaluating matrix elements of 
higher dimensional operators. 

The correlation function (|ip satisfies the standard dispersion relation 

ReHi 2(Q^) = - / ^^ds + subtraction. (11) 

TX J s + Q 

In QSR, Ren(Q^) in the left hand side of (jl^) is calculated by OPE at large as 
give in (^), while Imn(s) in the right hand side is parametrized by the nucleon pole 
and the phenomenological continuum. The pole part reads 

qlmlif\s) + lmnf\s) = 7r\%{q + M^)5{s - M^), (12) 

where Atv is defined as {0\ri\N) = X^u^p) with u{p) being the nucleon Dirac spinor. 
We assume that, when s > Sn, the hadronic continuum reduces to the same form 
with that obtained by an analytic continuation of OPE; 



Imn^°"*(s) = 7Te{s-SN) 



647r4 I TT V 12 V^ V J 327r2' tt 



(13) 

Imnr (3) = -ire{s - S^)^-^s (l + . (14) 

Substituting (^ and (0,|l^,|l^) into the dispersion relation (|IlD, one can gener- 
ate sum rules for the resonance parameters. We will write here the Borel transformed 
version of the sum rules in which the higher dimensional operators and the effect of 
the continuum can be relatively suppressed. (See Appendix A for useful formula.) 



2^ 



8 \ 12 ^ / 

6M2 , , al ,5 1 ,a,(M2)\ 

4i-Ai,M„e-"«/'"= = "^M^Erix) { 1 + l'^AMll\ _ ^. (i6) 

4 Y 2 7r y 72 

Here is the Borel mass, En = l — {1+X + -- - + ^)e"^ with x = Sn/M"^, 'Je is 
the Euler constant (0.5772- • •) and 
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ag = -47r'(gg), 6 = 47r'(-G"). (17) 

TT 

Note that we chose /i^ to be M ^ which is a typical scale of the system after the Borel 
transform. We call eq. (p!5D (eq. (|T6|)) as "even" ("odd") sum rule since it contains 
only even (odd) dimensional operators. As for the running coupling constant, we use 
a simplest one-loop form, 

°-<*^^' - 9.os(MVA-) - <^'* 

with = {0.l74:GeVy which is obtained to reproduce ^^(l) ~ 0.4 

Eq.(^) and Eq. ( [T6| ) will be used when we analyse the sum rules for g,,]^. One 
can derive formula for as a function of M in three ways; (i) the ratio of (p!5D and 
(0), (ii) the ratio of (^) and its logarithmic derivative with respect to and (iii) 
the ratio of (|16|) and its logarithmic derivative with respect to M^. We have made 
an extensive Borel stability analyses for these three cases. We found that the higher 
dimensional operator and the as corrections improve the Borel stability as well as the 
prediction for M^r. However, there is a wide range of which can reproduce the 
experimental nucleon mass in the Borel analyses. In later sections, we will use the 
finite energy sum rules (FESR) to fix S^- 

In Fig.l, shown are the Borel curves for three cases (i)-(iii) with = 
imiGeV^ obtained by the FESR with {qq) = -{0.2402GeV)^ and (^G^) = 
0.012GeV'^. The solid, dashed and dash-dotted curves correspond to cases (i), (ii) 
and (iii), respectively. The curve in case (iii) is the most stable of three cases and 
reproduces the experimental value. The others do not show good stability. Also, 
three curves are rather sensitive to the change of S^. 



III. SUM RULES FOR THE tt - N COUPLING CONSTANT 



In this section, we look at the two-point correlation function, 

KpiQ) = ^ J d'xe^'^-{0\TUx)U0)\7c\p = 0)), (19) 

where |vr°(p = 0)) is a neutral pion state with vanishing four momentum (the soft 
pion). Since we are working in the chiral limit, this soft pion is simultaneously on 
shell. 

Since we are taking the matrix element between the vacuum and the soft pion 
in (p!9D, only the pseudo-scalar operator qij^q survives in OPE for Tri{x)f]{0). The 
pseudo- vector operator q^^f/^^q vanishes since (0| 57^755 |7r°(p)) oc — > 0. Thus the 
relevant terms in OPE read as follows: 
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W{Q^) = Ildim3{Q^) + T^dim5{Q^) + T^dim7{Q^), (20) 

ndi,ri3(<5 ) = -«75 ^^^^ \ 2~) ^ ' 

^dim^{Q^) = 0, (22) 

In- 

n,.„7(Q2) = ^75T^(-G'2)(0M*75C^|vr). (23) 

/i^ dependence of as is again implicit in the above equations. 

The diagrams corresponding to (PT|), ( P^j ) and ( P3[ ) are Fig. 1(d), 1(e) and 1(f) 
respectively. The Wilson coefficients with 0(0;^) corrections in (|21|), ( [2^ ) and ( p3D 
turn out to be identical to (H), (P) and ([ToD respectively. This can be explicitly 
checked by carrying out OPE with the background field method in the fixed point 
gauge. An alternative way to see this is the plain wave method. As an illustration, 
let's consider dimension 3 operator (Fig. 1(d)) and expand eq.(|l]) in terms of dd and 

I J d^xe"'-''T{T]4x)r]p{0)) = Cs{l)apdd + Cps{il^)o.^di^d + ■■■ (24) 

where Cg^ps is the Wilson coefficient with s{ps) denoting scalar (pseudoscalar) . Note 
that ■ ■ ■ are vector, pseudovector and tensor operators. Note also that uu and ui'y^u 
do not arise for the loffe current. We sandwich eq.(p^) by free quark states to extract 
the Wilson coefficients. Applying the projections 1^/3 and (i75)a/3 to eq . (|2^) , one gets 



Cs = ^'^^p'l^ip^"^^ L, = y d'xe^'^-{p\T{r^{x)m)\p) (25) 
C,s = U,, = ^ / ciW-(p|T(r^(x)z75r](0))b) (26) 



{p\di-f,d\p) ' 4 

with \p) being a quark state with momentum p. Since massless QCD does not flip 
chirality in perturbation theory, {p\dd\p) ~ {p\di'~^^d\p) and Li ~ Lj^g except for 
trivial kinematical factors. Thus Cg = Cps is obtained even when as corrections are 
included. 

n'^((5^) satisfles the dispersion relation 

1 r Imll'^f s) 

ReWiQ'^) = - / ^ds + subtraction. (27) 

n J s + 

The hadronic imaginary part Imn'^(s) has a nucleon pole and the continuum. The 
pole part is parametrized by assuming the n — N — N vertex Cmt = gnNNi'y^T ■ nN 
as was done in . The continuum part is extracted from the analytic continuation 
of OPE. In total, 
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Imn"(s) = i^5TrX%g^NSis - Ml) 

-^757r^ (l + ^ v) ' S^md^ld\A. (28) 

where S-^ is the continuum threshold. Note that S'vr does not have to be the same 
with ^iv- 

Putting Eq.(^) and Eg. (pOD into (pT]) and making the Borel improvement, one 
obtains the following sum rule: 

"''^1(4) fl + 1^1 - i-UiA . (29) 




47r2 V 2 7r / 18 ' vr 

where is the pion decay constant (93 MeV), and {Q\di'-)^d\TT^) has been rewritten 
by the soft pion theorem 

(0|(ii75rf|7r°(p = 0)) = ^{dd). (30) 

One can get rid of the coefficient e^^/^^^ /X]^ in (^ by using the nucleon sum rules 
(p!5|) or (0). Thus one arrives at two different sum rules for (yf^^r: 



^£2(14) (i + + + li£o(|f4) + Tf (i - (i - hi.)^ 



which is obtained from the "even" sum rule (|T5|) for the nucleon, and 



(31) 



which is obtained from the "odd" sum rule ([Tq) for the nucleon. For M^v in eg . (|3^) , 
we just take the experimental number instead of reexpressing M^v by the Borel mass 
M through the nucleon sum rule. Even if one uses M7v(M) in eg. (|32D , the Borel curve 
for g°^^ is not affected so much as far as one adopts the odd nucleon sum rule (case (iii) 
in section II). Note here that, if one assumes 5*7^ = Sn, Eg. (|3^) gives g°^^ = M^/ f^,, 
namely the GT relation with qa = ^■ 

It is in order here to remark some difference of the present work from that of 
ref. 0,^. First of all, we have carried out OPE up to dimension 7 both for g^^fyf and 
the nucleon mass, while only the lowest dimension operator, namely qi'jbq, is taken 
into account for the g^N sum rules in ref. . Secondly, «s corrections to the Wilson 
coefficients are not taken into account in @J^. Thirdly, the continuum thresholds S'^r 
and Sjy are completely neglected in [^,|| for the g^^^ sum rules. 
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IV. BOREL ANALYSIS FOR THE n - N COUPLING CONSTANT 



A. The determination of the thresholds 3.,^ and Sn 

As we have mentioned in section II, it is difficult to determine Sn from the 
Borel sum rules of the nucleon. In this section, we will utilize the finite energy sum 
rules (FESR) to give a constraint on Sn as well as 5^. 

The FESR has a general form \T^\ , 



fSo fSo 

/ ds ImUoPEis) = dss"" ImUphenXs), (33) 
Jo Jo 

where lmIlopE{s) is an imaginary part obtained by the analytic continuation of 11 in 
OPE, Imllphen.i^) is the phenomenological imaginary part and 5*0 is the continuum 
threshold (either Sn or S^)- For the nucleon, by using Hi with n = 0, 1 and 112 with 
n = 0, one gets p 

64.'A=, = (l + ^ + 2.^(^G^>S„ + (l - 1^) . (34) 

64x'A^Af„ = -Sn^gg) (l + Sl + ^^'{gg){^G\ (35) 

QA.'XlMl = (l + f ^) ^ + .\'^G')Sl - ^-^Amr^SN. (36) 

V07r/4 TT 9 TT 

The renormalization point here is chosen as /i^ = Sn- 

By using the standard values of the condensates (gg)(lGeV^) = —(225 ± 
25MeVY and (^G^) = 0.012GeV\ we solved Eq.(|4D ~ Eq.(|6D numerically. Table 
1 shows the results for four different values of {qq). 





{-0.250GeVf 


{-OMOGeVf 


(-0.225Ge1/)3 


(-0.200Ge1/)3 


SN{GeV^) 


1.77 


1.60 


1.34 


0.887 


\N{GeV^) 


0.0267 


0.0235 


0.0187 


0.0118 


MN{GeV) 


0.997 


0.940 


0.845 


0.645 



Tablel: Sn, ^n, Mn obtained from Eq.(|3^) ~ (p6| ) with four different values of (qq). 
(qq) = — (0.2402GeV^)^ reproduces the nucleon mass. 



From this table, we choose Sn = 1-34 — 1.77 as physical range where the 
nucleon mass is reasonable reproduced within ±10% errors. Our Sn is smaller than 
that usually used in the literatures p-§,|T^. However, the as corrections are not 
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taken into account in these references. The effect of the correction to the spectral 
parameters can be exphcitly seen by expanding the solutions of Eq. (|^) up to linear 
in as] 

Xl = Xl(l- 26.2i-16n^qq)rH-G')) (1 - 3.08^), (37) 

Sn = S%(i- 21.2{-16n'{qq))-l{^G')) (1 - 3.26^), (38) 
\ vr / vr 

n / / o , , N -1 , ttq „o, \ , etc 

fU I 1 1 o «/' 1 «„^/;^„\N-o / * /^^\ * - " . -s 



Mat = M]; 1 - 18.6(-167r^(gg))-3(^G'0 (1 - 1.94—), (39) 

V TT / TT 

where Aq = 4(gg)^, S"^ = (6407r'^(gg)^)^ and = (—^Tc'^iqq))^, which are the 
solutions when the corrections and the gluon condensate are neglected. (|37D -(p9D 



show that the as corrections tend to reduce the observables by considerable amount 
particularly in Sn- 

Next, we estimate by taking the n = FESR of fl" 



Tvr. 



Xl,9,. = { (l H- 1 - ^{?G=> \ I ~{dd} I . (40) 




Since the FESR is rather sensitive to the structure of the continuum compared to 
the Borel sum rule and we do not know much about the detailed structure of the 
continuum for (|l^), we just limit ourselves to the n = sum rule (local duality 
relation) for safety. To roughly evaluate the range of S',r, we simply put g-^N = 13.4 
in (|40|) with A^r being determined in the nucleon FESR. The result is given in Table 
2 for three different values of the condensate: 





{qq)iGeV^) 


SNiGeV^) 


S^{GeV^) 


set 1 


-(0.250)3 


1.77 


1.98 


set 2 


-(0.240)3 


1.60 


1.85 


set 3 


-(0.225)3 


1.34 


1.62 



Table 2: SnjSt^ with three different values of {qq). 5*^ is obtained by substituting 

^^=13.4 into (|0D. 



From Table 2, one finds that St^ is always greater than Sn- This is consistent 
with the Borel sum rule g°j^ in ( p2D which tells us that Qt^n > M^/ f-n only when 
St^ > Sn- In subsections below, we will examine the Borel stability of Qt^n with the 
parameter sets obtained in Table 2. 
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B. Borel analysis for 5^*^" 

In Fig. 3, g^^'^ is shown as a function of M^. The sohd, dashed and dash-dotted 
curves correspond to set 1, set 2 and set 3 in Table 2, respectively. g'^§^ in Fig. 3 (a) 
includes the corrections to the Wilson coefficients, while they are neglected in 
Fig. 3(b) except for the gluon condensate. 

^e?jen j^g^g ^ sizablc variation and good Borel stability is not seen in Fig.3. 
By comparing Fig.3 (a) with Fig.3 (b), one finds that the corrections improve the 
Borel stability only slightly. The effect of the higher dimensional operator to the Borel 
curve is also small. In fact, + f^)}, which is a ratio of the 

dimension 3 term and the dimension 7 term in eg. (pQ]) , is about 4 % at ~ IGeV^. 
(Note that dimension 5 terms do not arise for the loffe current.) 

Although g'^^^ in eq.(PT]) is proportional to {qq), three curves in Fig.3, which 
correspond to different values of (gg), almost overlap with each other. This is because 
the change of (gg) is compensated by the changes of S'^^at. 

In Fig. 4, g'^§'^ for set 2 is shown with S.,^ being changed by ±10%. The solid, 
dashed and dash-dotted curves correspond to St, = 1.85 x 1.1, 1.85 and 1.85 x 0.9, 
respectively, f?^/^" increases as S'^ increases, which is consistent with the prediction 
of FESR in eq.(|g). 



C. Borel analysis for g°j^ 

Fig. 5 (a), (b) show ^f""^ as a function of M^. The solid, dashed, dash-dotted 
curves correspond to set 1, set 2, set 3, respectively, g^j^ in Fig. 5(a) includes as 
corrections to the Wilson coefficients, while they are neglected in Fig. 5(b) except for 
the gluon condensate. 

has apparently better Borel stability than g^)^"^, but still sizable varia- 
tion is seen. The corrections do not affect the Borel stability much, since the same 
as correction appears both in numerator and denominator in eg . (^) . 

g°j^ is rather sensitive to the change of the parameter sets , in particular the 
Stt/Sn- To see the effect of St^^n on g°j^ in more detail, we expand eq.(|32|) up to 

\ 1 

C2 b c2 b I ( 
^TT 9 '^N gj ) 

Sn)] , (41) 



,,,^,M^ (g^-§) f 2 
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where we have neglected small corrections and used the fact S*^^ ^ 6/9. The 
approximate formula eg .(PI) is in good agreement with the exact one eq. (|32D in 10 % 
for M2 > l.OGeVl 

When > Sn, independent term in (|41|) gives g^j^ = {S^/ SnY^Mn / f^) 
which is larger than M^/ f-w- The leading 1/M^ correction reduces g'^^ slightly. The 
experimental value for g^^^ = 13.4 is obtained when = 1.6GeV^ for the parameter 
set 3. 

In Fig. 6, g"^^ for set 2 is shown with 5*^ being changed by ±10%. The solid, 
dashed and dash-dotted curves correspond to = 1.85 x 1.1,1.85 and 1.85 x 0.9, 
respectively. g°j^ increases as S'^r increases, which is consistent with the prediction of 
FESR in eq. (|40|) and also with the approximate formula (pl). 



D. Comparison of g^^^"' and g°j^ 

As we have already mentioned, the Borel stability for g°'^j^ is better than g^)^"'. 
This is consistent with the fact that the Borel curve for M^r is most stable in "odd" 
sum rule (case (iii) in Fig.2). Also the absolute value of g"^^ is larger than f?^*^" 
for appropriate range of M^: e.g. g°j^ =11.3 — 12.8 versus g^)^"^ =9.02 — 9.10 at 
M2 = IGeVl 

In previous subsections, 5*7^ = l.85GeV^ has been used as a standard value to 
study the Borel stability. Alternative way to make the Borel analysis is to eliminate 
from eq.(|31|) and eq.(|3^) using eq.(|40D, and then to solve (7,rAf self-consistently for 
each M^. By this procedure, we found that there exists no solution satisfying eq.(^) 
and eq.(^OD simultaneously, while there exists solutions of eq.(]3^) and eq.(^OD which 



are given in Fig. 7. This result again confirms that g°]^ is better starting point to 
study the n — N coupling constant than f?^^". 



V. CONCLUSION 

In this article we have made extensive Borel and FESR analyses of gj^N by 
taking into account higher dimensional operators, as corrections and the continuum 
threshold. None of them has been considered in the previous analyses which led to 

g^r, = Mj,/U ii- 

What we have found are summarized as follows: 

(a) The higher dimensional operators up to dim. 7 and the as corrections play no 
crucial role for the Borel stability of (^^at. 

(b) (7°*^ is more appropriate for examining (yf^^v than g'^^'^-, since the former has better 
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Borel stability. This fact is also consistent with the fact that "odd" sum rule forM^y 
has a best stability. 

(c) Q-kN is most sensitive to the ratio 5"^/ Sjq, and both the FESR and Borel sum rules 
tell us that St,/Sn > 1 is a crucial ingredient to reproduce the experimental q-kn^ 

We also found that the Borel stability of 5'°'^, even if dim. 7 operator is taken 
into account, is not satisfactory enough to determine the pion-nucleon coupling con- 
stant precisely. For the nucleon sum rule, there have been some attempts to improve 
the Borel stability such as the modification of the loffe current [14| and the inclusion 
of the instantons |]TB|]. In particular, instantons improve the nucleon Borel sum rules 
considerably at low region, so it will be an interesting problem to study Qt^n with 
instanton contribution in the future. 
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^This point may have some relation to the Adler-Weisberger sum rule |16| which tells 
us that QA > ^ (or equivalently Qt^n > Mj^/f-,^) is obtained only when the continuum 
contribution in the tt — N channel is taken into account. 
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APPENDIX A: 



The Borel transform is defined as, 



with — * oo, n ^ oo, and ^ = being fixed. 

When B is apphed to the correlation function Il{Q^) = ^ J ^^^q2 , it leads to 

BU(Q^) = ^— / lmn(s)e"^ (A.2) 

This shows that the Borel transform tends to suppress the high energy contribution. 
Some useful formula are 

B{Q^)HogQ^ = {-l)''+^T{k + 1)(M2)^ (A.4) 

B(\ogQY = -21ogM2 + 27B, (A.7) 
SQ2(logg')' = 2M2(logM' - 7B + 1), (A.8) 
S(g2)2(iogg2)2 ^ (M2)2(_41ogM2 + 47s - 6). (A.9) 
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Figure Captions 

Fig.l 

OPE up to dimension 7 operators for the correlation of loffe current. Wavy lines 
denote gluon lines, broken lines denote the quark/gluon condensate. 
Fig. 2 

Mn (nucleon mass) as a function of the Borel mass squared M^. The solid, 
dashed, dash-dotted lines correspond to the cases (i), (ii) and (iii), respectively. 
{qq) = -{0.240GeV)^ {^G^) = 0M2GeV^ and Sn = l.QOGeV^ are used. 
Fig.3(a),(b) 

5'^'^" as a function of M^. The solid, dashed and dash-dotted lines correspond to set 
1, set 2 and set 3, respectively, corrections are taken into account in Fig. 3 (a), 
while they are neglected in Fig. 3(b) except for gluon condensate. 
Fig.4 

with {qq) = -(0.240Gel^)^ as a function of M^. The solid, dashed and dash- 
dotted lines correspond to S,, = 1.85 x 1.1, 1.85 and 1.85 x 0.9, respectively. 
Fig.5(a),(b) 

g°j^ as a function of JVP. The solid, dashed and dash-dotted lines correspond to set 1, 
set 2 and set 3, respectively, as corrections are taken into account in Fig. 5(a), while 
they are neglected in Fig. 5(b) except for gluon condensate. 
Fig.6 

5^°^ with (qq) = — (0.240GeV)^ as a function of M^. The sohd, dashed and dash- 
dotted lines correspond to 5^ = 1.85 x 1.1, 1.85 and 1.85 x 0.9, respectively. 
Fig.7 

(^ttTV as a function of M^. S*^ in eq.(^Dp with A at being determined in the nucleon 
FESR is used for S*^ in eq. (p2D . The solid, dashed and dash-dotted lines correspond 
to gf^ with (qq) = -{0.25GeVy, -{0.2A0GeVy and -{0.225GeVy, respectively. 
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